PLC-mediated PI(4,5)P 2 hydrolysis regulates activation and inactivation 
of TRPC6/7 channels 
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Transient receptor potential classical (or canonical) (TRPC)3, TRPC6, and TRPC7 are a subfamily of TRPC chan- 
nels activated by diacylglycerol (DAG) produced through the hydrolysis of phosphatidylinositol 4,5-bisphosphate 
(PI(4,5)P 2 ) by phospholipase C (PLC). PI(4,5)P 2 depletion by a heterologously expressed phosphatase inhibits 
TRPC3, TRPC6, and TRPC7 activity independently of DAG; however, the physiological role of PI (4,5) P 2 reduction 
on channel activity remains unclear. We used Forster resonance energy transfer (FRET) to measure PI(4,5)P 2 or 
DAG dynamics concurrently with TRPC6 or TRPC7 currents after agonist stimulation of receptors that couple to 
G q and thereby activate PLC. Measurements made at different levels of receptor activation revealed a correlation 
between the kinetics of PI(4,5)P 2 reduction and those of receptor-operated TRPC6 and TRPC7 current activation 
and inactivation. In contrast, DAG production correlated with channel activation but not inactivation; moreover, 
the time course of channel inactivation was unchanged in protein kinase C-insensitive mutants. These results sug- 
gest that inactivation of receptor-operated TRPC currents is primarily mediated by the dissociation of PI(4,5)P 2 . 
We determined the functional dissociation constant of PI (4,5) P 2 to TRPC channels using FRET of the PLCS Pleck- 
strin homology domain (PHd), which binds PI(4,5)P 2 , and used this constant to fit our experimental data to a 
model in which channel gating is controlled by PI(4,5)P 2 and DAG. This model predicted similar FRET dynamics 
of the PHd to measured FRET in either human embryonic kidney cells or smooth muscle cells, whereas a model 
lacking PI(4,5)P 2 regulation failed to reproduce the experimental data, confirming the inhibitory role of PI(4,5)P 2 
depletion on TRPC currents. Our model also explains various PLC-dependent characteristics of channel activity, 
including limitation of maximum open probability, shortening of the peak time, and the bell-shaped response of 
total current. In conclusion, our studies demonstrate a fundamental role for PI(4,5)P 2 in regulating TRPC6 and 
TRPC7 activity triggered by PLC-coupled receptor stimulation. 



INTRODUCTION 

Transient receptor potential classical/canonical (TRPC) 
3, C6, and C7 channels are the closest mammalian ho- 
mologues of the Drosophila melanogasterTRP channel and 
are expressed in various cell types, including smooth 
muscle and neurons (Hardie, 2003; Inoue et al., 2006; 
Venkatachalam and Montell, 2007) . These channels con- 
duct cations (Na + , Ca 2+ ) in response to stimulation of 
receptors coupled to phospholipase C (PLC), namely 
G q protein-coupled receptors and certain tyrosine ki- 
nase receptors (Ramsey et al., 2006). For that reason, 
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the currents mediated by these channels are often called 
receptor-operated cation currents (Inoue and Kuriyama, 
1993; Firth et al., 2007). TRPC3/6/7 channels are also 
activated by synthetic membrane-permeable diacylglyc- 
erol (DAG) analogues and are thus considered to be 
DAG-sensitive or activated channels (Hofmann et al., 
1999; Okada et al., 1999). In a physiological context, 
DAG is produced by the hydrolytic activity of PLC, 
which is located downstream of the receptors for neu- 
rotransmitters and hormones. Therefore, the receptor 
stimulation activates TRPC3/6/7 channels through the 
production of DAG to generate the receptor-operated 
TRPC currents (Beech et al., 2004; Panda et al., 2005; 
Hartmann et al., 2008) . 
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Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2, or 
PIP9) is the major substrate of PLC, and its hydrolysis 
produces DAG. PI(4,5)P2 is known to regulate numer- 
ous ion channels, modulating electrical signal outputs 
from metabotropic receptors in diverse physiological 
contexts (Gamper and Shapiro, 2007; Hilgemann, 2007; 
Logothetis et al., 2007). However, knowledge concern- 
ing the effect of PI(4,5)P2 on TRPC channels is still ac- 
cumulating (Jardin et al., 2008; Lemonnier et al., 2008; 
Monet et al., 2012). We have recently demonstrated 
using Danio rerio voltage-sensing phosphatase (DrVSP) 
that reduction or depletion of PI (4,5) Po inhibits the ac- 
tivity of DAG-sensitive TRPC3/6/7 channels, both in an 
exogenous expression system and in smooth muscle- 
derived cells (A7r5) (Imai et al., 2012; Itsuki et al., 
2012a). The DrVSP-mediated inhibition of TRPC3/6/7 
currents was detected even in currents evoked by a 
membrane-permeable DAG analogue (OAG) , which sug- 
gests that reduction in PI(4,5)P 2 can inhibit TRPC3/6/7 
channel opening regardless of the presence of DAG. 
Under G q protein-coupled receptor stimulation, PI(4,5)P2 
hydrolysis (breakdown) by the receptor-activated PLC 
largely contributes to the production of DAG. There- 
fore, such a complex relationship of PI(4,5)P2 and DAG 
suggests that TRPC3/6/7 channel activity may be regu- 
lated in a self-limiting manner. 

The effects of PLC-coupled receptor-driven channel 
regulation via enzymatic hydrolysis of PI (4,5) P 2 are not 
known. To improve our understanding, we simultane- 
ously measured PI(4,5)P 2 or DAG dynamics and recep- 
tor-operated TRPC currents evoked by carbachol (CCh; 
a muscarinic receptor agonist) or vasopressin (a vasocon- 
strictor). To this end, we measured the levels of PI(4,5)P2 
and DAG using a quantitative Forster resonance energy 
transfer (FRET)-based sensor alongside detection of 
TRPC6 and TRPC7 currents, which are more sensi- 
tive than TRPC3 currents to reduction in PI(4,5)P2 in 
human embryonic kidney (HEK)293 cells and smooth 
muscle-derived cells. We found that the temporal FRET 
dynamics of PI(4,5)P2 reduction closely correlate with 
the time course of activation and inactivation of receptor- 
operated TRPC6/7 channel currents. We also constructed 
a kinetic model of receptor-driven PI(4,5)P2-DAG signal- 
ing. This model, calibrated with DrVSP-derived func- 
tional dissociation constants for PI(4,5)P2 binding to 
TRPC3/6/7 channels, closely resembled our experi- 
mental results. Our study, combining experimental and 
computer simulation data, revealed the crucial role of 
receptor-stimulated PI(4,5)P2 hydrolysis in TRPC6/7 cur- 
rents. Part of the work presented here has appeared in 
abstract form (Itsuki et al., 2012b). 

MATERIALS AND METHODS 

Plasmids and cells 

The pcDNA3 expression vector encoding human TRPC6 (Gen- 
Bank accession no. NM_004621) was provided by T. Hofmann 



(Institut fur Pharmakologie und Toxikologie, Zurich, Switzer- 
land); pCI-neo expression vectors encoding mouse TRPC3 (Gen- 
Bank accession no. NM_019510) and TRPC7 (GenBank accession 
no. NM_012035) were provided by Y. Mori (Kyoto University, 
Kyoto, Japan). Single amino acid mutation in TRPC6 and TRPC7 
was generated using the QuikChange Site-Directed Mutagenesis 
kit (Agilent Technologies) according to the manufacturer's in- 
structions. To generate bright FRET pairs, super-enhanced YFP or 
CFP isolated from a RI10A FRET sensor (provided by M. Matsuda, 
Kyoto University, Kyoto, Japan) were modified to give A207K mu- 
tants as a monomeric form (CFP mse or YFP mse ) (Zacharias et al., 
2002) . These modified fluorophores were fused to the N-terminal 
side of the PLCS Pleckstrin homology domain (PHd; provided by 
K. Jalink, The Netherlands Cancer Institute, Amsterdam, Nether- 
lands) to construct PI(4,5)P 2 sensor molecules consisting of CF- 
P mse -PHd or YFP rase -PHd. For DAG detection, CFP mse was fused to 
the C-terminal side of PKCe (provided by M. Schaefer, Leipzig 
University, Leipzig, Germany), yielding PKCs-CFP mse . To be an 
energy acceptor of membrane-bound PKCE-CFP mse , YFP mse was 
attached to the C-terminal side of the GAP-43 myristoyl domain 
(Invitrogen) through an octaglycine (G8) linker (Myr-YFP mse ). 
PI(4,5)P2 and DAG sensor cDNA were each incorporated into an 
IRES-reporter region-excluded pIRES2 expression vector (Invit- 
rogen). A pEF-BOS expression vector encoding human musca- 
rinic type 1 receptor (MjR) was provided by T. Haga (Gakushuin 
University, Tokyojapan). Human phosphatidylinositol-4-phosphate- 
5-kinase (PIP5K; |3 isoform) in pcDNA3.1 vector (Invitrogen) was 
provided by S. Kita and T. Iwamoto (Fukuoka University, Fukuoka, 
Japan). All PGR products were sequenced entirely. 

HEK293 cells (obtained from the ATCC) were maintained in 
Dulbecco's modified Eagle's medium (Invitrogen) supplemented 
with 10% FBS (Gibco) and antibiotics (penicillin and streptomy- 
cin; Gibco) at 37°C (5% C0 2 ) . For transfection, cells were seeded 
on poly-L-lysine-coated glass coverslips (Matsunami) in 35-mm 
culture dishes and transfected with a mixture of plasmid vector- 
incorporated DNAs using the SuperFect transfection reagent 
(QIAGEN). For FRET-based PI(4,5)P 2 detection, HEK293 cells 
were cotransfected with 1 pg each of plasmids encoding TRPC3, 
6, or 7 together with MjR (or without it, in endogenous musca- 
rinic receptor stimulation) and 0.3 pg each of plasmids encoding 
CFP mse -PHd and YFP mse -PHd. For DAG detection, 0.3 pg each of 
plasmids encoding Myr-YFP mse and PKCe-CFP mse were cotrans- 
fected instead of the PI(4,5)P2 sensor plasmids. For detection 
of local PI(4,5)P 2 around the TRPC7 channel, the sequence en- 
coding the donor protein (CFP mse ) was inserted before the stop 
codon of TRPC7. In this case, equal amounts (1 pg) of donor, ac- 
ceptor (YFP mse -PHd) , and MjR plasmids were used for transfec- 
tion. Measurements on transfected cells were made within 24—72 h 
after transfection. 

A7r5 cells, the cell line derived from rat thoracic aortic smooth 
muscle (Brandt et al., 1976), were obtained from the ATCC, 
maintained in medium identical to that used for HEK293 cells, 
and passaged every 5-7 d. The transfection protocol was essen- 
tially the same as the one used with HEK293 cells. A7r5 cells trans- 
fected with CFP mse -PHd and YFP mse -PHd were reseeded on poly-L- 
lysine-coated glass coverslips and incubated at 37°C (5% C0 2 ) 
for at least 15 min before use. Cells were always used within 2 h 
of reseeding. 

Solutions and drugs 

The standard external solution contained (iriM): 140 NaCl, 5 KC1, 

1 CaCl 2 , 1.2 MgCl 2 , 10 HEPES, and 10 glucose (pH 7.4, adjusted 
with Tris base; 300 mOsm, adjusted with glucose). The pipette 
solution contained (mM): 120 CsOH, 120 aspartate, 20 CsCl, 

2 MgCl 2 , 5 EGTA, 1.5 CaCl 2 , 10 HEPES, 2 ATP-Na 2 , 0.1 GTP, and 
10 glucose (pH 7.2, adjusted with Tris base; 290-295 mOsm, adjusted 
with glucose). CCh (Sigma-Aldrich) was diluted in the standard 
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external solution from its stock concentration (100 mM in H 2 0). 
To confirm monovalent cationic currents, NMDG solution (150 mM 
JV-methyl-D-glucamine chloride, 10 mM HEPES, and 1 mM CaCl 2 , 
with pH 7.4 adjusted with HC1) was applied at the end of each 
stimulus. RHC80267 (EMD Millipore) was dissolved in DMSO 
(Wako Chemicals USA) . Stock solutions of Arginine 8 vasopressin 
(AVP; 100 uM; MP Biomedicals) and nifedipine (10 mM; EMD 
Millipore) were dissolved in H z O and DMSO, respectively. AVP 
and nifedipine were freshly prepared in the standard external so- 
lution to final concentrations of 1 and 5 pM, respectively, before 
applying to A7r5 cells. During experiments, HEK293 and A7r5 
cells were continuously perfused with external solution and gravity- 
fed at a flow rate of 0.25 ml/ min. The perfusion was turned on and 
off using electromagnetic solenoid microvalves (The Lee Co.). 

Simultaneous measurements of TRPC currents and FRET 

Electrophysiology. The whole-cell patch-clamp technique was used 
for current detection. Patch electrodes with a resistance of 4—6 MCI 
(when filled with internal solution) were made from 1 .5-mm boro- 
silicate glass capillaries (Sutter Instrument). Series resistance 
errors were compensated >60%. Voltage generation and current 
signal acquisition were accomplished vising a patch-clamp amplifier 
(AxoPatch 200B; Axon Instruments) with an A/D D/A converter 
(Digidata 1200; Axon Instruments). Sampled data were low-pass 
filtered and digitized at 1 kHz using pClamp 9.0 (Axon Instru- 
ments) and analyzed using custom-written software (MATLAB; 
MathWorks) . The currents were recorded at a holding potential 
of —50 mV. For activation of DrVSP, depolarizing step pulses 
(from 20 to 180 mV, 500-ms duration) were delivered every 20 s. 
The ratio of the currents before and after DrVSP activation was 
used to quantitate DrVSP-mediated inhibition, "r (I) ." Before its 
calculation, the leak that was defined by the current in NMDG- 
containing solution was subtracted. All experiments were performed 
at room temperature (22-25°C). 

FRET detection. Fluorescence from voltage-clamped cells was 
detected using a microscope (60 x 0.9 N.A. objective; TE300 Eclipse; 
Nikon) equipped with a two-channel simultaneous beam-splitter 
(Dual-View2; Photometries) and a high sensitivity EMCCD cam- 
era (EvolveS 12; Photometries) . Excitation light filtered at 427/10 
and 504/12 nm was alternately introduced via an optical fiber 
from a lamp house equipped with a high speed excitation wave- 
length selector (75 W xenon lamp; OSP-EXA; Olympus) . Epifluo- 
rescence from the cells was prefiltered using a multiband dichroic 
mirror (449-483 and 530-569 nm) contained in the microscope, 
and then further separated in the beam-splitter (at 505 nm) and 
filtered at 464/23 nm (detection of the donor fluorescence) or 
542/27 nm (detection of the acceptor fluorescence). Optical fil- 
ters were obtained from Semrock, except the splitter (Chroma 
Technology Corp.). The duration of camera exposure was 100 ms 
and occurred within 150-ms periods of illumination at each exci- 
tation wavelength. Images were captured with an EM gain of 300 
and then digitized as 5 1 2 x 5 1 2 pixels by 1 6-bit arrays in the micro- 
scope software (Micro-manager v. 1.4) . The image pixel resolution 
was ~0.26 urn. Averaged intensities from the whole-cell region 
(typically 20 x 20 to 40 x 40 square pixels) were analyzed to calcu- 
late FRET using a custom-written MATLAB program. The electro- 
physiology and FRET measurements were synchronized using 
brief triggers from the A/D D/A converter linked to the excita- 
tion light shutter. All of the data in this paper were recorded from 
the first application of any of the agonists. 

Calculation of FRET. Fluorescence signal output obtained from a 
given sample is denoted by the descriptor iipo> where X and Y are 
the fluorescence filter settings for the emission and excitation light, 
respectively. The emission filter for donor fluorescence (464 nm) 



is denoted as "464," and that for acceptor fluorescence (542 nm) 
as "542." The excitation filter for donor excitation (427 nm) is 
denoted as "D," and that for acceptor excitation (504 nm) as "A." 
Background intensity, captured using the corresponding filter 
setting with nontransfected cells, was subtracted from specimen 
fluorescence signals. Finally, the FRET ratio (FR) was calculated 
according to the "3-cube" method (Erickson et al., 2001): 

*R = (W) -SD1- F 464(D) )/EA- (F 542(A) - i?D2 • F 542(D) ) , ( 1 ) 

where EDI = ^(dj/Fb^d), MJ2 = F 542(A) /,F 464(D) , and RA = 
^*542(D)/^542(A)- Constants of EDI, ED2, and EA were predeter- 
mined using measurements from single cells expressing only 
donor- (CFP mse ) and acceptor- (YFP mse ) tagged molecules, respec- 
tively. From the FR values, we can compute the effective FRET ef- 
ficiency (£ E ff): 

E„ = (FR - 1) • [ £ YFP mse (427) / £ CFP mse (427)], (2 ) 

where eYFP mse and eCFP mse are the molar extinction coefficients 
for the FRET cube excitation filters obtained using the 427-nm 
excitation band-pass filter. We determined the ratio in brackets to 
be 0.11, based on maximal extinction coefficients forYFP mse and 
CFP mse (Morietal., 2011). 

Establishment of a relationship between PI(4,5)P 2 
concentration and FRET 

Here, we calculated the relationship between FRET and PI(4,5)P 2 . 
To this end, we considered two different cases: ( 1 ) FRET between 
TRPC7 channel-anchored CFP mse and YFP mse -PHd, where YFP mse - 
PHd can transit between a soluble cytoplasmic state and a PI (4,5) P 2 - 
bound membrane state; and (2) FRET between CFP mse -PHd and 
YFP mse -PHd, where both probes transition between a soluble and 
membrane-bound state. 
Case 1 is described by: 



FRmo* 1 l K <l(PHd) (3) 

[P/(4,5)P 2 ] 

where F is the fraction of PHd bound to PI(4,5)P 2 , K diPHJ) is the 
dissociation constant of PHd bound to PI(4,5)P 2 (reported as 
2.0 pM; Lemmon etal., 1996; Hirose etal., 1999), andFfi„„is the 
maximum FR at an infinitely high concentration of PI (4,5) P 2 that 
induces all the fluorophore-fused PHd probes to bind to the 
plasma membrane. FE max is a purely theoretical value, because 
physiological PI(4,5)P 2 levels (5-40 pM) are insufficient to local- 
ize all PHd proteins to the plasma membrane (Bunce et al., 1993; 
McLaughlin and Murray, 2005) . Nevertheless, although FR max is 
impossible to demonstrate, overexpression of PIP5K can greatly 
increase the cellular levels of PI(4,5)P 2 by approximately two- to 
threefold (Winks etal., 2005). We observed that cells overexpress- 
ing PIP5K demonstrated ~1.2 times higher FR than control cells 
with resting PI(4,5)P 2 levels (Fig. 5 A). We thus calculated the 
FE„ mx value by multiplying the resting FE by this correction factor 
of 1.2. Eq. 3 was used to simulate FR dynamics in cells express- 
ing TRPC7-CFP mse and YFP mse -PHd, as shown in Eq. 25 (Table 3) 
and Fig. 7 D. Conversion of the FRET between CFP mse -PHd and 
YFP mse -PHd to PI(4,5)P 2 concentration (i.e., case 2) is described 
in the Results. 

Estimation of expressed fluorophore-tagged PHd proteins 

We determined that the average concentration of fluorophore- 
tagged PHd in single cells was 1.6 pM. This value was calculated 
based on the intensity from fluorescein (Sigma-Aldrich) as follows. 
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Fluorescein was dissolved in 50 mM borate buffer (pH 9.1). The 
fluorescence intensities in small droplets of this fluorescein solution 
were measured, under the same conditions as the cells overex- 
pressing YFP mse -PHd. The respective intensities were standardized 
by the quantum yields for fluorescein (0.92) and YFP mse (0.57), 
and then the average [YFP mse -PHd] in living cell was determined 
to be 0.8 pM. Equal amounts of CFP mse -PHd and YFP mse -PHd plas- 
mids were transfected, so total amounts of CFP mse -PHd and 
YFP mse -PHd proteins could be extrapolated from this value. 

Modeling and fitting statistics 

Kinetic models of PI(4,5)P 2 hydrolysis pathway by PLC were for- 
mulated as ordinary differential equations, and the concentra- 
tions of PI(4,5)P 2 and DAG derived from these were incorporated 
into the channel-operation models. Simulations were performed 
in Excel (Microsoft) using the forward Euler method with a time 
step of 0.05 s. Individual steps were translated into differential 
equations based on the proposed kinetic scheme. Fitting the mod- 
els to the experimental data of TRPC6/7 currents was performed 
by a Generalized Reduced Gradient algorithm of the Solver func- 
tion in Excel. Details of model formulations and model fitting to 
the experimental data are described in the Results. 

The errors between the experimental and the back-calculated 
i<Rby model fitting to the currents (Figs. 6 and 7) were evaluated 
by the standard deviations of residual as follows: 

SD(o-) = ^i; =0 (Norm.Fi? - Norm.***). 7(» - 1), 

where Norm.i^R and Norm.i i K wi denote the normalized experi- 
mental and back-calculated FRat the i points, n indicates the total 
number of time points (20/ s) . 

Online supplemental material 

Plotting the currents versus FR^,,/ FR,^,^ relationship (Fig. SI). 
DrVSP-mediated inhibition at the single-channel level (Fig. S2). 
Diffused PI(4,5)P 2 was incorporated in the self-limiting regula- 
tion by PI(4,5)P 2 -DAG signaling (SPD) model-based simulation 



(Fig. S3). Incomplete matching of the DG model to the experi- 
mental FR (Fig. S4). The quick recovery of FR under the AVP 
stimulation (Fig. S5). Fitting the SPD model to TRPC7 current 
demonstrated less matching to FR dynamics and vice versa 
(Fig. S6) . Effect of an inactive ATP analogue in the patch pipette 
on the receptor-operated TRPC7 currents (Fig. S7). The online 
supplemental material is available at http://www.jgp.org/cgi/ 
content/full/jgp.201311033/DCl. 

RESULTS 

Simultaneous measurement of PI(4,5)P 2 and TRPC currents 

The PHd of PLCS binds both PI(4,5)P 2 and inositol 
1,4,5-trisphosphate (jp 3 ) (Hirose etal., 1999). In the rest- 
ing state, however, most of the fluorophore-tagged PHds 
are located at the plasma membrane, making it possible 
to measure PI(4,5)P 2 at the plasma membrane (van der 
Wal et al., 2001; Jensen et al., 2009; Yudin et al., 2011) 
(Fig. 1 A) . Using this approach, we first demonstrated the 
simultaneous measurements of PI(4,5)P 2 and receptor- 
operated TRPC6 current. The FRET pairs, which con- 
sisted of donor (CFP mse ) or acceptor (YFP mse ) fused to the 
PHd, were coexpressed with TRPC6 channel and IV^R 
in mammalian HEK293 cells (Fig. 1 B). Fig. 1 C shows a 
typical example of the simultaneous measurement of a 
TRPC6 current and PI(4,5)P 2 levels in a HEK293 cell 
after stimulation with 10 pM CCh. Soon after CCh ap- 
plication, an inward TRPC6 current (Fig. 1 C, top) and 
FRET reduction (middle) were observed concurrentiy. 

To evaluate the kinetics of the current activation and 
inactivation, we measured the time required for recep- 
tor-operated TRPC6 current to increase from 10 to 90% 
of its peak amplitude and then to decay from 90 to 50% 
(Fig. 1 C, top) . We quantitatively examined the real-time 
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Figure 1 . Simultaneous measurement of receptor-operated 
TRPC6 currents and PI(4,5)P 2 detected by 3-cube FRET. 

(A) Images of PI (4,5) P 2 sensor expressed in HEK293 cells 
using a 3-cube filter. Phase-contrast image (top left), YFP 
channel (F M2 (a)> top right), CFP channel (F iMiB y, bottom 
left), and FRET channel (i"54 2 (D); bottom right) are shown. 

(B) Diagram of molecules transfected into HEK293 cells. 
TRPC6, PI(4,5)P 2 sensor (CFP mse -PHd [blue box] and 
YFP mse -PHd [yellow box]), and IV^R were expressed. The 
excitation wavelengths, 427 and 504 nm, were alternately 
illuminated. (C) Typical example of CCh-induced TRPC6 
currents (top) and the corresponding FRET changes (mid- 
dle). The FR (middle; circles) calculated by 3-cube meth- 
ods can yield near absolute FRET efficiency (E E¥¥ ; right 
axis). Measured parameters were the duration of A10-90% 
and A90-50% of the peak current, the kinetics of FRET 
decay (tFR), and minimum FR (FR^n). The decline of 
FRET (FR) was fitted with a single-exponential decay (red 
solid curve): FR = FR^ + FR A x exp(-t / tFR) . The bot- 
tom panel shows the changes in the fluorescence intensi- 
ties (a.u.) that passed through the respective filter setting. 
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alteration in PI (4,5) P? level using the 3-cube FRET mea- 
surement (described in Materials and methods). FRET 
changes over time between CFP mse -PHd and YFP mse -PHd 
(Fig. 1 C, middle) were calculated from the respective 
fluorescence intensities (Fig. 1 C, bottom). The donor 
fluorescence (Fig. 1 C, open circles; i^tD)) and reso- 
nance fluorescence upon donor excitation (black circles; 
^342(D>) displayed inverted changes, whereas acceptor 



fluorescence (filled triangles; F^^a)) stayed largely con- 
stant (Fig. 1 C, bottom). These configurations ensured 
the FRET reduction and low level of quenching of sen- 
sor proteins during the recordings. Kinetics of FR re- 
duction (tFR) and the minimum amount of FR (-FRmin) 
under the receptor stimulation were obtained by fitting 
to the FRET data with a single-exponential decay func- 
tion (Fig. 1 C, middle). 
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Figure 2. Correlation between the TRPC6/7 currents and the decay of PI(4,5)P 2 . (A and B) Example traces of TRPC6 currents (top) 
and FRET of PI(4,5)P 2 sensor (CFP mse -PHd and YFP mse -PHd) (bottom) upon stimulation with 100 uM CCh of either endogenous 
(A) or overexpressed IVfR (B). (C and D) Same as in A and B, but in cells expressing TRPC7. (E) Summary of currents and FRET 
changes. TRPC6/7 current increase (A10-90%) and decay (A90-50%), kinetics of FRET reduction (tFK), and degree of FRET reduc- 
tion (FRmin) were accelerated in a CCh concentration and IVfR expression-dependent manner. The data depicted by the stripe and 
the white bars show without channel expression (transfected only PI(4,5)P 2 sensor) and endogenous receptor simulation, respectively. 
Numbers in parentheses indicate the number of cells measured, here and throughout. (F and G) Time courses of the A10-90% and 
A90-50% of receptor-operated currents were plotted against the simultaneously measured iFR (F, TRPC6; G, TRPC7). These data were 
obtained from various concentrations of CCh or level of IVfR expression. The slope with a linear fit highlights a relationship between 
the time courses and tFR 
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PI(4,5)P 2 dynamics at different levels of receptor 
stimulation and TRPC6/7 currents 

We then explored the effect of PI(4,5)P 2 reduction on 
TRPC6 or TRPC7 channel currents at varying levels of 
receptor stimulation. Stimulation of endogenous mus- 
carinic receptors in HEK293 cells with 100 pM CCh 
evoked prolonged TRPC6 current and a small amount 
of reduction in FR (Fig. 2 A) . The time for TRPC6 cur- 
rent activation (A10-90%) and inactivation (A90-50%) 
were 32 ± 7 s and 39.2 ± 6.7 s (n = 6) , respectively. In ac- 
cordance with the current time course, jFRv/as length- 
ened to 55.3 ± 4.5 s, and FR^^^ stayed near the resting 
level. Reduction in i<Rwas only 18 ± 6%. Previous stud- 
ies have shown that depletion of PI(4,5)P 2 can be in- 
duced by overexpression of M : R (Xie et al., 2011; Dickson 
et al., 2013). Consistent with these studies, overexpress- 
ing MjR and treating the cells with a high concentration 
of CCh (100 pM) greatly accelerated the TRPC6 current 
and FRET reduction (Fig. 2 B) (A10-90% = 1.6 ± 0.3 s; 
A90-50% = 6.2 ± 0.9 s; tFR = 4.5 ± 0.9 s), and enhanced 



the FR reduction with a near zero FRET efficiency (FR nin = 
1.15 ± 0.1 and £ EFF = 0.016 + 0.01; n = 11). Intriguingly, 
the time to reach the peak current was remarkably 
shortened compared with endogenous muscarinic re- 
ceptor stimulation at the same concentration of CCh 
(100 pM; +MiR = 2.8 ± 0.7 s; endo = 63 ± 8 s). Similar 
tendencies in the FR reduction and the peak time were 
observed when TRPC7 was expressed instead of TRPC6, 
except for the shorter activation and inactivation time 
for TRPC7 current (Fig. 2, C and D) . In addition, TRPC7's 
tFR was slightly delayed and FRm^ was also slightly at- 
tenuated compared with that of TRPC6. (Data from the 
different strength of receptor stimulation was summa- 
rized in Fig. 2 E.) 

To elucidate the functionality of PI(4,5)P 2 hydrolysis, 
we focused on the kinetic relationships between TRPC6/7 
currents and FRET reduction. For that purpose, the log- 
log plots for the activation or the inactivation of TRPC6/7 
currents and the values of tFR were made using data ob- 
tained at varying levels of agonist stimulation, with and 
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Figure 3. Incompatible correlation between receptor-operated TRPC6/7 currents and DAG production. (A) Principle of DAG detec- 
tion by PKC probe FRET. Increments in FRET caused by the translocation of PKCe-CFP mse to the plasma membrane were detected by 
a coexpressed membrane-resident acceptor protein (Mry-YFP mse ) in HEK293 cells. (B) Whole-cell TRPC6 currents (top) and FRET 
changes caused by DAG increments, "FiT 1 ^" (bottom, green circles), recorded from endogenous muscarinic receptor stimulation with 
100 pM CCh. The rise of FRET was fitted to the exponential equation: FR A " g = FR d * s mlIX - AFi? dag x exp(-t / rFR^A (green solid curve). 
(C) Traces of currents and FI& g in M^-overexpressing cells with 100 pM CCh. (Left) TRPC6 currents. (Right) TRPC7 currents. Pro- 
longed DAG production was observed. Purple zones indicate inconsistencies between current inactivation and DAG production. The 
inset in the TRPC7 panel shows FFr^ changes over a longer time scale (300 s). (D) Summary of FFr^ g levels at the respective current 
points observed in the robust receptor stimulation (+MiR and 100 pM). The black and gray bars denote expression of TRPC6 and 
TRPC7 channels, respectively. (E) Time courses of initial phase of current increase (A10-30%) were plotted against kinetics of DAG 
production (tF^). 
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without MiR overexpression (Fig. 2, F and G) . These 
log-log plots showed a clear correlation between tFR 
and the current time courses (A10-90% [left panels] and 
A90-50% [right panels]) for both TRPC6 and TRPC7 
currents. The linear relationship that appeared in 
a plot of the relation of A90-50% to tFR in TRPC7- 
expressing cells showed significantly steeper slopes 
(slope = 1.33) than that in TRPC6-expressing cells (slope = 
0.76) (Fig. 2, F and G) . This steepness may reflect 
higher TRPC7 sensitivity to reduction in PI (4,5) P? than 
TRPC3 or TRPC6 (Imai et al., 2012). In addition to the 
FRET decay, the extent of the reduction or depletion 
of PI(4,5)P2 levels (-FRmm) also showed a similar ten- 
dency to the time course of current activation or inacti- 
vation of each channel (Fig. SI). These plots, however, 
were slightly more scattered than those of tFR. Such a 
scattering was probably caused by cell-to-cell variability in 
the released IP3 in response to the hydrolysis of PI(4,5)Ps 
(Irvine and Schell, 2001). The simultaneous detection of 
PI(4,5)P9 and TRPC currents demonstrated that the ac- 
tivation and inactivation time courses related to both 
the kinetics and the extent of PI (4,5) P2 reduction. 

Simultaneous detection of DAG and receptor-operated 
TRPC channel current 

TRPC3/6/7 channels are DAG-sensitive ion channels, 
but the manner in which DAG dynamics correlate with 
these TRP channels' activity remains largely unknown. 
To address this question, DAG production was con- 
currently monitored with receptor-operated TRPC6/7 
channel currents. The detection of DAG dynamics re- 
lies on membrane translocation of DAG-activated PKC 
in response to increasing DAG levels at the plasma 
membrane (Violin et al., 2003) (Fig. 3 A) . We used Ca 2+ - 
insensitive PKCe as a fluorescence donor molecule to 
exclude Ca 2+ -dependent translocation of PKC (Sinnecker 
and Schaefer, 2004) . Upon stimulation with CCh through 
the endogenous muscarinic receptors, TRPC6 current 
and DAG production were initiated almost simultane- 
ously (Fig. 3 B) . This parallel response is consistent with 
TRPC6 being a DAG-sensitive channel. Furthermore, 
during the inactivation of the TRPC6 channel current, 
DAG level also declined. This synchronicity indicated 
that when the strength of receptor stimulation was weak, 
the production of DAG levels seemed to be a critical fac- 
tor to the current appearances. 

Contrary to the synchronicity of the TRPC6 current 
and DAG dynamics in the weak receptor stimulation, 
when the cells overexpressed MiR, the simultaneous 
measurement was revealed to be inconsistent. The acti- 
vation of TRPC6 channel current paralleled DAG pro- 
duction, whereas the inactivation of the channel did 
not; there was no decline in DAG production (Fig. 3 C, 
left, purple zone). This lack of temporal consistency 
between the current decay and DAG levels was even 
more prominent in the TRPC7 channel, which exhibited 



inactivation while DAG levels were still increasing (Fig. 3 C, 
right, and summarized in D) . This inconsistency may be 
explained by the idea that PKC-mediated phosphoryla- 
tion inhibits channel opening, which has been pro- 
posed for TRPC3 (Soboloff et al., 2007). 

However, for TRPC6 channel, there was no significant 
difference in the current inactivation (A90-50%) be- 
tween cells overexpressing PKCe (TRPC6, MjR and DAG 
sensor-expressing, 7.5 ± 2.5 s; n = 5) and control (TRPC6 
and MiR, 6.2 + 0.9 s; n = 11). In addition, a previous 
report has shown that PKC8 exerts a negative feedback 
effect via phosphorylation of Ser448 of TRPC6 (Bousquet 
et al., 2010). Using the phosphorylation-insensitive mu- 
tant TRPC6 S448A and its corresponding mutant TRPC7 S394A , 
we tested whether the PKC-mediated phosphorylation 
of these residues is involved in the current decay. Under 
the robust receptor stimulation, unexpectedly, we did 
not detect any clear difference in the A90-50% time 
(100 uM CCh: TRPC6 wt , 11.4 ± 2.4 s; TRPC6 S448A , 12.7 + 
2.7 s; TRPC7 wt , 1.8 ± 1.1 s; TRPC7 S394A , 1.9 ± 0.9 s; n=8). 
Furthermore, the basal phosphorylation site of TRPC6 
has been identified at Ser 814 (Bousquet et al., 2011), 
which may play a role in the channel function. We also 
tested the mutant, TRPC6 S814A , but it failed to show any 
significant differences, including A90-50% (10.1 ± 3.5 s; 
n = 7) , which is consistent with the previous report 
(Bousquet et al., 2011). Our data do not exclude the 
possibility that PKC-mediated phosphorylation inhibits 
channel opening, because of the variety of cell and 
measurement conditions. Nevertheless, these results 
suggest that the DAG production level and PKC phos- 
phorylation-mediated channel modulation may contrib- 
ute less to the inactivation of TRPC6/7 channels at the 
robust stimulation. 

A simulation model was built for further understand- 
ing of dual regulation by PI(4,5)Po and DAG, as described 
in the latter section. For that purpose, we presented the 
relation between the current increase and the produc- 
tion of DAG. Plotting the early phase of the current 
increase (A10-30%) versus DAG production kinetics 
(tFR 1 ^) exhibited a clear correlation, with a smaller slope 
for the TRPC7 compared with the TRPC6-expressing cells 
(Fig. 3 E). This result suggests that the TRPC7 channel 
is highly sensitive to the increment of DAG. We thus 
set the DAG sensitiveness in the initial parameter 
as TRPC7 > TRPC6 in the model simulation (Table 4, 
rows 21 and 22). 

Functional dissociation constants of PI(4,5)P 2 binding 
to TRPC3/6/7 channels 

These results have shown the substantial importance 
of the dissociation of PI(4,5)P9 to the inactivation of 
TRPC6/7 channel currents, but the affinity of PI(4,5)P 2 
to these channels is not yet known. To obtain this pa- 
rameter, DrVSP, which functions as a membrane-resident 
voltage-controllable phosphoinositides phosphatase, was 
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used to reduce intrinsic PI(4,5)P2 (Okamura et al., 
2009). Our previous report demonstrated that reduc- 
tion of PI(4,5)P2 by activation of DrVSP led to concomi- 
tant inhibition of TRPC3/6/7 currents (Imai etal., 2012). 
We simultaneously measured the voltage-dependent 
stepwise controls of FRET between CFP mse -PHd and 
YFP mse -PHd, and the DrVSP-mediated inhibition. To evoke 
the currents, the DAG lipase inhibitor, RHC80267, was 
used. This compound is suitable to produce stabilized 
inward currents, mainly caused by elevating the resting 
level of DAG (Albert et al., 2005) . The step-pulse proto- 
col, from 20 to 180 mV with a duration of 500 ms, en- 
ables the current inhibition and FRET reduction to be 
observed simultaneously (Fig. 4 A) . 

By plotting the current inhibition "r(7)" and FRET re- 
duction "r(FR)" against the depolarizing pulses that ac- 
tivate DrVSP, the various sensitivities of the TRPC3, 
TRPC6, and TRPC7 channels were quantified (Fig. 4 B) . 
The inhibition of the TRPC3 current (Fig. 4 B, left, circles) 
was relatively insensitive, compared with the reduction in 
FRET (triangles). In contrast, the inhibition of TRPC7 
was highly sensitive to the reduction in FRET (Fig. 4 B, 
right) . TRPC6 exhibited a similar level of gradual changes 
in r(I) and r(FR) (Fig. 4 B, middle) . This result confirms 
our previous observations showing differential sensitivi- 
ties of TRPC3/6/7 channels to PI(4,5)P 2 reduction, 
with an order of TRPC7 > C6 > C3 (Imai et al., 2012) . 

However, the inhibition of TRPC3 and TRPC6 currents 
was not complete, and the FRET reduction by DrVSP 



activation was also insufficient to reach a zero FRET 
level (FR = 1.46 ± 0.07; n = 14 at 120 mV). We thus spec- 
ulated that this partial current inhibition may simply be 
because of the incomplete depletion of PI(4,5)P 2 . This 
idea raises a challenging question: how is TRPC3/6/7 
inhibition by DrVSP activation observed at the single- 
channel level? To answer this question, DrVSP was acti- 
vated during receptor stimulation in the cell-attached 
patch mode. Robust depolarization led to a brief, but 
almost complete, inhibition of the TRPC6 channel dur- 
ing the bursting activity (Fig. S2) . This observation en- 
abled us to use a standard ligand-binding isotherm for 
PI(4,5)P 2 -TRPC3/6/7 channel binding. The calcula- 
tion of PI(4,5)P2 concentration from the reduction in 
FR was done by a boundary function (described in 
Materials and methods) . Because the FRET between 
CFP mse -PHd and YFP mse -PHd proteins is cancelled by 
detachment of either the donor or acceptor fluoro- 
phore-fused PHd proteins from the membrane, FR 
reduction could be approximated as a cooperative 
square law of the membrane-bound fraction of PHd 
(Eq. 3) as follows: 

-^,F 2 =fl/l+ feg> Tf, (4) 
FR max { { [P/(4,5)P 2 ]JJ 

where K d (p H<t) is the dissociation constant of PI(4,5)P2 
binding to the PHd. The left side of Eq. 4 is always posi- 
tive, and solving PI(4,5)P2 yields: 
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Figure 4. TRPC current inhibition and PI(4,5)P 2 
reduction in response to the protocol for measur- 
ing the voltage dependence of DrVSP activation. 
(A) TRPC6, CFP mse -PHd, YFP mse -PHd (PI(4,5)P 2 
sensor), and voltage-sensing phosphatase (DrVSP) 
were coexpressed in HEK293 cells. Gradual current 
inhibition and reduction in PI(4,5)P 2 caused by the 
step-pulse protocol (left; from 20 to 180 mV; dura- 
tion of 500 ms; repeated every 25 s) . A DAG lipase 
inhibitor (RHC80267; 100 uM) was applied to in- 
duce the currents (gray horizontal bar). The ratio 
of current inhibition, r(I), and FRET reduction, 
r(FR), upon the depolarization pulses was used to 
quantify the channel activity and PI(4,5)P 2 changes 
after DrVSP activation (right). (B) The voltage de- 
pendence of current inhibition (left axis; circles) 
and FR reduction (right axis; triangles) after DrVSP 
activation in cells expressing TRPC3 (left), TRPC6 
(middle) , and TRPC7 (right) channels. 
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Figure 5. Functional dissociation constants of PI(4,5)P 2 bind- 
ing to TRPC3/6/7 channels. (A) Comparison of FR or Eeff in the 
resting condition between cells expressing TRPC6 channel and 
CFP mse -PHd and YFP mse -PHd (control), and those overexpressing 
PIP5K (+PIP5K) . The FR of cells overexpressing PIP5K increased 
on average by ^1. 2-fold compared with control cells. *, P < 0.05; 
unpaired <test. (B) Steady-state plots for estimating the functional 
Ka of PI(4,5)P 2 binding to TRPC3/6/7 channels. Horizontal axis 
indicates the estimated PI(4,5)P 2 concentration based on the con- 
version from FR to PI(4,5)P 2 , according to Eq. 5. 



where FR^.^ is the maximum FR value at an infinite con- 
centration of PI(4,5)P 2 . According to the resting FR of 
PIP5K-overexpressing cells, that value was estimated as 
1.2-fold higher than the resting FR of the control cells 
(Fig. 5 A). After solving Eq. 5, r(I) versus the estimated 
PI(4,5)P 2 concentration (Est, PI(4,5)P 2 ) plots were fitted 
using the Hill equation (Hill, 1910) (Fig. 5 B). Assum- 
ing a dissociation constant of PI(4,5)P2 — PHd = 2 uM 
(Hirose et al., 1999), the functional dissociation con- 
stants of PI (4,5) P 2 binding to TRPC3, TRPC6, and TRPC7 
channels were estimated at 1, 2, and 5 uM, respectively. 
These factors were incorporated into K 3 as the initial 
parameter for simulation in channel regulation by 
PI(4,5)P 2 (see below). 

Modeling TRPC channel activity coupled to 
PI(4,5)P 2 -DAG signaling 

Having determined the functional constants of PI (4,5) P 2 
binding to TRPC channels, we attempted to simulate 
channel activity and compare the results with the experi- 
mental data. Our model consists of three components. 
Part 1 covers the minimal PI(4,5)P 2 -to-DAG reaction, 
including the process of PI(4,5)P 2 recovery (Table 1, 
Eqs. 6-17). Part 2 is the calculation of the open proba- 
bility (P 0 ) and the resultant current based on the dy- 
namics of DAG and PI(4,5)P 2 concentrations (Table 2, 
Eqs. 18-20). Part 3 is a back-calculation of normalized 
FR of PI(4,5)P 2 sensor according to the concentrations 
of PI(4,5)P 2 and IP 3 (Table 3; Eqs. 21-26) . Details of the 
respective components are described in full below. 
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the experimental data and elucidate the functional role 
of the reduction in PI(4,5)P9 caused by agonist-induced 
receptor-operated currents. Thus, we focused on the ac- 
tivation of PLC by G q protein-coupled receptors as fol- 
lows. The initial step is hydrolysis of PI (4,5) P? by PLC, 
which generates DAG and IP3 (Fig. 6 A, k\, and Table 1 , 
Eqs. 6-8). Activation of PLC starts at the first time step 
(0.05 s) after time zero. In the factor of PLC activity, we 
also included the adjusted factors of receptor desensiti- 
zation (Table 1 , Eq. 9) and the time-dependent accelera- 
tion of PLC activity by receptor stimulation (Eq. 10) . The 
latter factor has been demonstrated as Ca 2+ -dependent 
positive feedback in PLC activity (Horowitz et al., 2005) . 
The next step is DAG phosphorylation by DAG kinase to 
generate phosphatidic acid (PA) (Fig. 6 A, ftjj, and Table 1 , 
Eq. 11). The third step is production of PI(4)P, a pre- 
cursor of PI (4,5) P2, from PA, skipping the intermediate 



products of CDP-DAG and phosphatidylinositol (Fig. 6 A, 
hn, and Table 1, Eq. 12). For recovery of PI(4,5)Ps, we 
considered two pathways: (1) resynthesis of PI(4,5)Po 
from PI(4)P by PIP5K (Fig. 6 A, h» and Table 1, Eq. 13) 
and (2) diffusion of PI (4,5) P2 from global to local chan- 
nel area (Eqs. 14 and 15). The idea and necessity of the 
diffusion pathway are described in the Discussion and 
the legend of Fig. S4. The amount of scavenged PI(4,5)P2 
arising from the binding of PHd proteins was solved by 
Eq. 16 (Table 1). The total available amount of PI(4,5)P 2 
to transfer to Part 2 (channel gating) was calculated by 
Eq. 17 (Table 1). 

(Part 2) Channel gating. Two channel gating models were 
constructed. The first, a simple channel gating model, 
only consists of DAG binding and unbinding for channel 
opening and closing, and it lacks the inhibitory regulation 



TABLE 1 

Equations for minimal PI( 4,5)Pr-DAG reaction scheme 



Equation No. 
d[PI(4,5)P 2 ] M =d[Pl(4,5)P 2 ] tai -d[PI(4,5)P 2 ] irai . V De fw - dt 

d[DAG] = d[DAG] + d[Pl(4,5)P 2 ] K -£ i(1) - De t([) • dt-d[DAG] \ - dt 

d[lP 3 ] = d[IP 3 ] + d[Pl(4,5)P 2 ] Mt -A,;,, - De t(t) - dt-d[lP 3 ] • dt 

De f(t) =l-(Rd_f ex P (-x rd )/t) + Sd_f • exp(-x sd /t) 
Rd _ f : fraction of rapid desensitization 
Sd _ f : fraction of slow desensitization 
T rd : time constant for the rapid desensitization 
T sd : time constant for the slow desensitization 

*,(,) =K ■ B_PLC pg " 

B_PLC : basal PLC activity 

pg : power parameter of G q PCR receptor stimulation 

d[PA] = d[PA] + d[DAG] • h A - dt -d[PA] • dt 

d[PI(4)P] = d[PI(4)P] + d[PA] • k a - dt-d[PI(4)P] ■ £„■ dt 

d[Pl(4,5)P 2 ] rec _ d =d[PI(4)P] -K ■ dt-d[Pl(4,5)P 2 ] ie _ d . A, (l) . De f(t) . dt 

d [ PI (4,5)P 2 ] gl<iM =d[Pl(4,5)P 2 ] g|oM -d[PI(4,5)P 2 ] g|oM . W De f(t) . dt 
/fj g : rate constant for global PLC activity 

4 PI ( 4 ' 5 ) P 2 ]„ n ^ 

• (l-errf(spot/sprt(4 • dcoef • dt)))-d [PI(4,5)P 2 ] • k- i(l) - De f(t) - dt 
errf : error function 

dcoef : diffusion coefficient of PI (4,5) P 2 
spot : distance between local and global domain 

d [ PI (4,5)P 2 l„ engcd PHd =(d[PHd] re , -d[PHd], p3 ) • (l/(l + X dJHd /(d[PI(4,5)P 2 ] toi 
+ d[Pl(4,5)P 2 ] r _ d+ d[PI(4,5)P 2 ]„ n ))) 

d[PI(4,5)P 2 ] Bt =d[Pl(4,5)P 2 ] M+ d[Pl(4,5)P 2 ]^+d[Pl(4,5)P 2 ] difi ^ ta -d[Pl(4,5)P 2 ] 
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TABLE 2 

Equations for channel gating 



Comment to the equation 
Open probability in DG model: 

Open probability in SPD model: 



TRPC6 / C7 currents: 



Equation 

P 0 =(jf,-A' 2 /d[DAG] 2 +K 2 /d[DAG] + l)~' 

IK, -K. 2 / d [DAG] 2 + K 2 / d [DAG] + (if, -K 2 / d[DAG] 2 ) ■ (K 3 / d[Pl(4, 5)P 2 ] M ) 
° +K 2 - (x 3 /d[Pl(4,5)P 2 ] to[ )/d[DAG] + (X3/d[Pl(4,5)P 2 ] M ) + l 
^and^are dissociation constants for DAG to TRPC channel. 
K 3 is dissociation constant for PI(4,5)P 2 to TRPC channel. 

I = N-g-P 0 - (V h -V m ) 

N : number of channels 
g : single channel conductance 
V h : holding potential (-50 mV) 
: reversal potential (0 mV) 



No. 

18 



1!) 



1>0 



by PI (4,5) reduction (termed the "DG model" in Fig. 6 A, 
bottom left) . The second is a realistic model that takes 
into account the inhibitory effect caused by a reduction 
of PI(4,5)P2 that direcdy induces transition of the chan- 
nel to an inactive state from any closed or open state. 
This model is referred to as the "SPD model" (Fig. 6 A, 
bottom right). Open probabilities (P Q ) expressed as a 
function of agonist concentration, calculated according 
to the DG and SPD models, are described in Eqs. 18 and 
19, respectively (Table 2). Finally, whole-cell currents 
were calculated using Eq. 20. 

(Part 3) Back-calculation of FR from the simulated PI(4,5)P 2 
concentration. As described in the previous section, the 
observed FR changes were converted to PI(4,5)P2 con- 
centrations, according to Eq. 5. Here, we attempted to 
do the opposite, using these equations to explore FRET 
dynamics. Receptor stimulation produces DAG and IP3 
after PLC-mediated hydrolysis of PI(4,5)P2- IP3 is highly 



diffusible in the cytoplasm and binds to PHd proteins 
with higher affinity (0.1 uM) than PI(4,5)P 2 (Hirose et al., 
1999) . The factor of FRET reduction caused by the bind- 
ing of IP3 to PI(4,5)P9 sensor was incorporated into the 
back-calculation of FR dynamics (Table 3, Eqs. 21-24). 
Although it was difficult to calculate the absolute effi- 
ciency of FRET, because of various expression levels of 
fluorophore-fused PHd, we could estimate the normal- 
ized FR (Norm ..fR) changes by the resting level of FR The 
alternations in Norm.FR, caused by CFP mse -fused TRPC7 
channel versus YFP mse -PHd and CFP mse -PHd versus YFP mse - 
PHd, were solved by Eqs. 25 and 26 (Table 3) , respectively. 

Insufficient matching of the DG model 

First, we examined whether a model without PI (4,5) P? 
regulation (DG model) was able to reproduce the exper- 
imental data of the simultaneous measurements. Fitting 
to the TRPC6 currents was accomplished by minimizing 
the sum of the squared errors with 19 free parameters, 



TABLE 3 

Equations for the back-calculated FR 



Comment to the equation 



Equation 



No. 



Fraction of IP 3 bound PHd (PHd): 

Concentration of 

membrane bound PHd: 
Concentration of cytosolic PHd 

upon PI(4,5)P 2 reduction: 

Concentration of IP3 bound PHd: 

Norm.fR between TRPC7-CFP m!e and 
YFP m , P -PHd: 



d(PHd D , 3b ) = l-d[PHd] Ip3 /[PHd] tot 

d [ pHd L m =[ pHd L- d [ pHd L 

d[ PHd L, = (l "I / (A' drai / d[PI(4,5)P 2 ] m + l))-[PHd] K 
d [ PHd ]iP 3 =(l/(^d.iP3 / d [l p 3 ] + l))- d [ p Hd] td 



Norm.fR = d (PHd IP3b ) ■ ( (l / (l + A' iPHi / d[PI (4, 5)P 2 ] w )) • [l / (l + K irm / [PI (4, 5) P 2 ^ ) J 
j^Pl(4,5)P 2 ] t : Pl(4,5)P 2 concentration at time zero 



21 

22 
23 
24 
25 



Norm.fR between CFP„ 
andYFP m , P -PHd: 



,-PHd 



Norm.fR = d (PHd IP3b ) • (l/ (l + A" dPHd / d [PI (4, 5) P 2 ; 



1/ 1 + Kd. 



,/[ p l(4,5) 



J)" 



20 
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the initial values of which are listed in Table 4. The fi- 
delity of the model was assessed by a similarity of FR 
between the experimentally measured FR and simu- 
lated FR which was obtained from the back-calculation 
of the resultant PI(4,5)P2 concentrations by fitting to 
the currents, according to the equations described in 
the section on model Part 3 and Table 3. 

When the currents were induced by the weak recep- 
tor stimulation through endogenous muscarinic recep- 
tors, the experimental FR was substantially compatible 
with the calculated FR from the simulated PI(4,5)P 2 
changes (Fig. S4 A; SD = 0.12). In contrast, fitting the 
currents in MiR-overexpressing cells to the DG model 
highly deviated from the experimental FR changes 
(Fig. S4 B; SD = 0.63). Furthermore, the DAG produc- 
tion was quite more transient than that observed in the 
PKC-based FRET dynamics (Fig. S4 B, bottom panel, 
dashed line). Therefore our fitting examination indi- 
cates that the DG model may be useful for mimicking 



the delayed receptor-operated currents, but it is not to- 
tally suitable for the rapid case. As we demonstrated, 
when the parameter for PLC activity (fa) was set to an 
accelerated kinetics (fa = 0.7) , a marked TRPC6 current 
inactivation emerged only in the SPD model but not in 
the DG model (Fig. 6 B, top). 

Fitting the SPD model to the experimental data 

The SPD model was then tested. The same basic fitting 
strategy as in the DG model was used. After fitting the 
SPD model to the experimental receptor-operated cur- 
rent data, the computed PI(4,5)Pa data were compared 
with the experimental FR dynamics. In this case, in ad- 
dition to the 19 parameters, the dissociation constant of 
PI (4,5) P? binding was incorporated at K 3 (Table 4, row 
23). The fitting of the SPD model to the whole-cell 
TRPC6 currents did indeed show overlapping of the 
dynamics of the experimental Norm.i 7 /? and the back- 
calculated Norm.i<R, with a smaller SD value than in the 



TABLE 4 



Parameters and initial conditions for the fitting 


Row 


Parameters (unit) 


Setting 


Initial value of TRPC6 endo/C6 + M^R/ 
A7r5/C7 + MjR 


Source or comments 


1 


PI(4,5)P 2 (uM) at resting 


Free 


20/20/20/20 


Bunce et al., 1993; McLaughlin and Murray, 
2005 


2 


/ji_PLC (s" 1 ) 


Free 


0.04/1/1/1 


Rational to PI(4,5)P 2 reduction with our data 


3 


tDAG kinase (s _1 ) 


Free 


0.03/0.03/0.08/0.03 


Rational to DAG changes with our data 


4 


fei_PA to PIP reactions (s _1 ) 


Free 


0.01/0.01/0.01/0.01 


Appropriate for PI(4,5)P 2 synthesis 


5 


/ji V _PIP5K (s~ ] ) 


Free 


0.1/0.1/0.25/0.1 


Appropriate for PI(4,5)P 2 resynthesis 


6 


&v_IP3 phosphatase (s _1 ) 


Free 


0.5/0.5/0.5/0.5 


Appropriate for IP 3 hydrolysis 


7 


Trd (S) 


Free 


5/2/1/4 


Appropriate for Receptor desensitization 


8 


T,d (s) 


Free 


50/10/5/10 


Same as above 


9 


Rd_f (no unit) 


Free 


0.5/0.5/0.5/0.5 


Same as above 


10 


Sd_f (no unit) 


Free 


0.5/0.5/0.5/0.5 


Same as above 


11 


spot (distance global to local; um) 


Free 


4/4/3/4 


5 times more than the diffusion coefficient of 
PI(4,5)P, 


12 


dcoefof PI(4,5)P 2 (pm 2 /s) 


Free 


0.8/0.8/0.8/0.8 


Golebiewska et al., 2008 


13 


Ratio of local k / global k 


Free 


1/4/2/5 


Approximate from the uneven FRET reduction 
(Fig. S3) 


14 


PI4P (uM) 


Fixed 


10/10/10/10 


Brown etal., 2008 


15 


Activation delay (no unit) 


Free 


0.01/0.001/0.001/0.001 


Appropriate for receptor activation 


16 


Activation power (no unit) 


Free 


0.5/0.3/0.7/0.3 


Same as above 


17 


Vrev (mV) 


Fixed 


0/0/0/0 




18 


Vhold (mV) 


Fixed 


-50/ -50/ -50/ -50 




19 


No. of channels 


Free 


100-7,000 


Appropriate for the current density 


20 


Channel conductance (pS) 


Fixed 


35/35/35/70 


Hofmann et al., 1999; Lemonnier et al., 2008 


21 


Ki (Ki for DAG1; uM) 


Free 


60/60/35/10 


Effective OAG concentrations are 10 to 100 pM 
in Hofmann et al., 1999; Okada et al., 1999; Imai 
etal., 2012 


22 


K 2 (Ai for DAG2; uM) 


Free 


30/30/10/10 


Same as above 


2.°. 


K, (AiforPI(4,5)P 2 ; uM) 


Free 


2/2/5/5 


This paper, Fig. 4 D 


24 


Expressed PHd (pM) 


Fixed 


1.6/1.6/1.6/1.6 


This paper, Materials and methods 


25 


_KdPI(4,5)P 2 ofPHd (pM) 


Fixed 


2.0/2.0/2.0/2.0 


Hirose et al., 1999 


26 


Ki IP 3 of PHd (uM) 


Fixed 


0.1/0.1/0.1/0.1 


Hirose et al., 1999 
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DG model (SD = 0.04; n = 4) . Furthermore, similar over- 
lapping was clearly evident in TRPC6 currents in MjR- 
overexpressing cells (SD = 0.05; n = 4; Fig. 7 B, middle). 
The fitting expressed a prolonged existence of DAG, 
persisting well beyond the current decay (Fig. 7 B, bot- 
tom, solid dashed line) . Such prolonged existence of 
DAG was seen experimentally after robust receptor 
stimulation (Fig. 3 B) and in recent work by Falken- 
burger et al. (2013). The correlation between the ex- 
perimental and modeled data supports the description 
of TRPC6 currents and PI(4,5)P2 dynamics in the SPD 
model (fitted parameters are summarized in Table 5). 

The SPD model was also examined for TRPC6/7-like 
currents in aortic smooth muscle-derived A7r5 cells 
(Brueggemann et al., 2006; Maruyama et al., 2006) . The 
PI(4,5)P 2 sensor proteins (CFP mse -PHd and YFP mse -PHd) 
were coexpressed in A7r5 cells without exogenous ex- 
pression of channels or receptors (Fig. 7 C, inset) . By 
applying AVP, TRPC6/7-like currents and FRET reduc- 
tion were observed similarly to the observations after 
stimulation of HEK293 cells with CCh (Fig. 7 C, top and 
middle). A noticeable difference was that PI(4,5)P2 re- 
covered quickly (5-20 s after the application of AVP). 



Such rapid FRET recovery is expected with the rapid 
desensitization of vasopressin receptors, and was repro- 
duced by coexpressing vasopressin type 1A receptors 
with TRPC6 or TRPC7 channel in HEK293 cells (Fig. S5) . 
When fitting these AVP-evoked TRPC6/7-like currents, 
the simulation parameters for rapid desensitization 
(Table 4, rows 7 and 8) also overlapped the Norm.i^K 
dynamics (Fig. 7 C, middle, blue line; SD = 0.12) . These 
results indicate that the SPD model is useful even in the 
context of physiological cells. 

Contrary to these positive results, the back-calculated 
FRET demonstrated less similarity to the experimental 
FRET achieved by fitting to TRPC7 currents observed in 
HEK293 cells and vice versa (Fig. S6; SD = 0.18; n= 4). 
In addition, TRPC7 currents often demonstrated a pla- 
teau or biphasic response when MiR was overexpressed, 
but the FRET did not clearly show such irregular dynam- 
ics. We reasoned that the FRET detection by fluorophore- 
fused PHd was too slow to respond to the PI(4,5)P2 
dynamics compared with the time course of TRPC7 cur- 
rent, recorded by the electrophysiological method. To 
improve this issue, we redesigned the FRET pairs to de- 
tect PI(4,5)P2 changes in the local vicinity of the TRPC7 
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Figure 7. SPD model fitting to experimentally observed TRPC6/7 currents and AVP-evoked TRPC6/7-like currents in A7r5 cells. 
(A) Fitting of the receptor-operated TRPC6 current simultaneously measured with PI(4,5)P 2 by FRET with the SPD model (top; light 
blue trace) . The experimental data were obtained from the low strength of receptor stimulation carried by CCh application (100 pM) 
to endogenous muscarinic receptor in HEK293 cells. Back-calculated FR of PI(4,5)P 2 concentrations at the respective points (bottom; 
solid) was normalized against the time zero FR (Norm.i<R; middle; light blue), and that was overlaid on the experimental Norm.FR 
(middle; circles). The bottom panel displays the fitting resultant changes in PI(4,5)P2 (solid), DAG (dashed), and PA (thin dashed). 
The similarity of Norm.i'T? was assessed by SD as described in Materials and methods. (B) Fitting of experimental TRPC6 current data 
from fy^R-overexpressing cells. The initial parameters are identical to those in A. (C) Fitting of experimental TRPC6/7-like currents 
recorded from A7r5 aorta-derived smooth muscle cells. The currents were evoked by 1 pM AVP. The inset in the top panel shows A7r5 
cells expressing the PI(4,5)P 2 sensor. (D) Fitting of TRPC7 currents. Local PI(4,5)P 2 dynamics were detected using a FRET donor linked 
to the TRPC7 channel at the end of its C-terminal domain (top; inset) . Transient increments in FRET were detected during the plateau 
or biphasic response (middle; arrow) . 
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channels. The donor fluorophore (CFP mse ) directly 
fused to the channel was coexpressed with YFP mse -PHd 
(Fig. 7 D, inset) . Intriguingly, this FRET pair showed 
a transient increment of FRET during the plateau or 
biphasic current (Fig. 7 D, red arrow). By fitting the 
TRPC7 currents with the corresponding "down-up- 
down" PI(4,5)P 2 dynamics (Fig. 7 D, bottom), we achieved 
an improved matching of Norm.FR (SD = 0.14) as well. 
This transient up-regulation of PI(4,5)P2 reflects a 
quick replenishment by diffusible PI(4,5)P9 in the SPD 
model (Fig. S3) . These results show that corresponding 
Norm.FR responses in various cell and channel settings 
strongly support the fidelity of the SPD model for re- 
producing simultaneous events of receptor-operated 
TRPC currents and FR changes. 

DISCUSSION 

Because PI(4,5)P9 hydrolysis by PLC is the critical event 
in receptor-operated TRPC3/6/7 currents, the parallel 
observation of PI(4,5)P2 or DAG with channel activity is 
fundamental to advancing our understanding of chan- 
nel regulation. To address this point, we performed 
concerted quantitative FRET measurements of sensitized 



fluorescence emission (3-cube method) and electro- 
physiological measurements on patch-clamped cells. 
We accomplished this by using CFP mse - or YFP mse -fused 
PHd proteins as a FRET sensor of PI(4,5)P 2 , or CFP mse - 
fused PKCe and YFP mse -fused myristoyl membrane- 
attached peptides to monitor DAG. CCh or AVP was 
used to induce TRPC6/7 currents, which were mea- 
sured to monitor the synchronicity of receptor-operated 
currents and changes in either PI(4,5)P9 or DAG levels. 
The experimental results demonstrated a correlation 
between the kinetics of PI(4,5)P2 reduction and the acti- 
vation or the inactivation of receptor-operated TRPC6/7 
currents. Fitting the experimental data to the DG or 
SPD model revealed that the idea of self-limiting regula- 
tion by PI(4,5)P2 and DAG is critical for reproducing 
receptor-operated TRPC6/7 currents. 

Reproducing the currents and FRET reduction 
in the SPD model 

Our results show that reduction in PI(4,5)P2 by PLG 
mediated hydrolysis appears to be fundamental to the 
inactivation of TRPC6/7 channels. This contribution is 
even more apparent after robust receptor stimulation, 
either by high concentrations of receptor agonist or by 



TABLE 5 



Resultant parameters by fitting to TRPC currents with SPD model 


Cell 




HEK293 


HEK293 


A7r5 


HEK293 


Transfected plasmids 


TRPC6/CFPmse-PHd/ 


TRPC6/M,R/CFPmse- 


CFPmse-PHd/YFPmse- 


TRPC7-CFPmse/M 1 R/ 






YFPmse-PHd 


PHd/YFPmse-PHd 


PHd 


YFPmse-PHd 


Receptor agonist (cone. uM) 


CCh (100) 


CCh (100) 


Arg- Vasopressin (1) 


CCh (100) 


The number of fitted cells 


4 


4 


4 


4 


Row 


Parameters (unit) 










1 


PI(4,5)P 2 (uM) at resting 


22.5 ± 0.4 


2.3.0 ± 0.4 


22.2 ±0.2 


14.6 ±0.1 


2 


k_I>LC (s _1 ) 


0.05 + 0.02 


1.2 ±0.20 


1.3 ±0.42 


0.93 ± 0.09 


3 


Mi_DAG kinase (s -1 ) 


0.04 ±0.01 


0.03 ± 0.02 


0.07 ± 0.0.3 


0.02 ±0.01 


4 


Mii_PA to PIP reactions, (s _1 ) 


0.008 ± 0.002 


0.006 ± 0.002 


0.004 ±0.001 


0.006 ± 0.0001 


5 


£ iv _PIP5K (s _1 ) 


0.050 + 0.004 


0.07 ±0.04 


0.17 ±0.083 


0.04 ± 0.03 


6 


*»_IPs phosphatase (s~ ] ) 


0.7 ±0.4 


0.2 ±0.1 


0.88 ± 0.093 


0.84 ±0.14 


7 


Trd (S) 


12.0 ±3.8 


1.4 ±0.6 


1.4 ± 0.2 


4.4 ±0.4 


8 


T,d (s) 


24 ± 11 


17 ±5.2 


6 ±2.4 


.38 ± 15 


9 


Rd_f (no unit) 


0.65 ± 0.08 


0.55 ±0.03 


0.60 ±0.12 


0.46 ± 0.05 


10 


Sd_f (no unit) 


0.35 ± 0.05 


0.44 ± 0.02 


0.43 ±0.124 


0.53 ± 0.06 


11 


spot (distance global to local; um) 0.64 ± 0.09 


2.78 ± 0.79 


2.8.3 ± 0.98 


1.90 ±0.29 


12 


dcoef ofPI(4,5)P 2 (um 2 /s) 


0.36 ± 0.08 


0.8 ±0.47 


0.5 ±0.11 


0.29 ± 0.06 


13 


Ratio of local h\ / global \ 


1.1 ±0.28 


4.0 ± 1.6 


2.8 ± 0.98 


6.0 ± 1.3 


14 


Activation delay factor 


0.0.3 ± 0.024 


0.01 ±0.002 


0.003 ± 0.002 


0.03 ± 0.025 


15 


Activation power factor 


0.6 ± 0.33 


0.32 ±0.07 


0.7 ±0.17 


0.3 ±0.12 


16 


The number of channels 


1,644 ±277 


5,808 ± 959 


317 ± 121 


2,384 ± 1,896 


17 


A'i (^forDAGl; uM) 


.39 ± 6.3 


26 ±2.4 


23 ± 7.6 


7.0 ± 0.3 


18 


K 2 (£ d forDAG2; uM) 


29 ± 5.1 


22 ± 2.0 


8 ± 3.4 


6.3 ± 0.25 


19 


K, (£dforPI(4,5)P 2 ;uM) 


1.9 ±0.6 


3.0 ± 0.66 


9 ±1.5 


9.7 ±0.14 




SD 


0.04 


0.05 


0.12 


0.14 




Sro.of 


82.7% 


80.9% 


70.1% 


56.2% 



Parameters are presented as mean + SEM. (Typical results are displayed in Fig. 7.) 
Percentage of the squared residuals lower than the value of 0.01. 
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overexpression of functional receptors. Here, we simu- 
late in a virtual setting the inhibitory role played by a 
reduction in PI(4,5)P2 levels in the SPD model in the 
context of different levels of PLC activity. In the aver- 
aged fitted data for TRPC6 currents, PLC activity (k) 
was gradually decreased or increased, as shown in Fig. 8 A. 
The simulated TRPC6 currents and FRET by the PHd 
clearly demonstrated PLC-dependent dynamics in re- 
gard to the amplitudes or kinetics in the currents (Fig. 
8 A, top) and the kinetics or the extent of FRET reduc- 
tion (bottom). We then depicted log-log plots of the 
time course of current growth (A10-90%) or decay 
(A90-50%) versus tFR (Fig. 8 B, red circles). For com- 
parison, we also simulated this plot with the DC model 
(Fig. 8 B, black circles) as well as TRPC7 channels. The 
log-log plots of A10-90% to tFR showed a linear rela- 
tionship in both SPD and DC models of TRPC6/7 cur- 
rents. These plots closely resemble the experimental 
data, which are shown as gray lines in Fig. 8 B (solid and 
dashed line for TRPC6 and TRPC7 channels, respec- 
tively) . In contrast, the linear relationship between cur- 
rent decay (A90-50%) and jFR was reproduced only 
when modeled with the SPD model (Fig. 8 B, right, red 
symbols) , but not with the DC model. 

The simulation result indicates that when the strength 
of receptor stimulation PI(4,5)P2 is low or the reduction 
of PI (4,5) P, is delayed (tFR >50 s), the DC model may 
be useful to reproduce such delayed receptor-operated 
TRPC currents. However, the situation alters as receptor 
stimulation is strengthened or the reduction of PI (4,5) Po 
is accelerated, and the DC model deviates from the 



experimental data. Furthermore, the SPD model-based 
relationship in A90-50% TRPC7 currents to tFR resulted 
in a slight increase in the steepness to TRPC6 (Fig. 8 B) . 
This result is similar to the TRPC7 current to iFR rela- 
tionship observed experimentally (Fig. 8 B, dotted 
line) . In contrast, the SPD model-based TRPC6 relation- 
ship resulted in an even steeper slope compared with 
the experimentally observed relationship (Fig. 8 B, solid 
line). The difference implies further unknown charac- 
teristics that underlie the regulation by PI(4,5)P2 and 
DAG in receptor-operated TRPC6/7 currents. Recently, 
several reports have addressed an alternation of PI (4,5) P2 
regulation by intracellular Ca' + in KCNQ2/3 (Kosenko 
et al., 2012; Falkenburger et al., 2013) and TRPV6 chan- 
nels (Cao et al., 2013). According to Kosenko et al. 
(2012), calmodulin binding to KCNQ2 channels stabi- 
lized PI(4,5)P2 channel binding. The involvement of 
Ca 2+ regulation to PI(4,5)P2 or DAG regulations requires 
further investigation. Nevertheless, because of the con- 
sistency of the SPD model simulation with the experi- 
mental results, we concluded that the inhibitory effect 
of reductions in PI (4,5) P2 would accelerate the channel 
inactivation with the increased receptor stimulation. 

Proposed function of the self-limiting regulation 
in receptor-operated TRPC6/7 currents 

Based on the model simulation, we explored several 
PLC-dependent features of the self-limiting regulation 
by PI(4,5)P2 and DAG. The first characteristic is the 
maximum P a (P omix ) . As shown in the left plot in Fig. 8 C, 
in both models, P om ax values gradually increased as PLC 
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Figure 8. Proposed contribution of PI(4,5)P2 
reduction in receptor-operated TRPC6/7 cur- 
rents. (A) Representative SPD model simula- 
tions for receptor-operated TRPC6 currents 
(top) and FRET by PI(4,5)P 2 sensor (bottom) 
from various PLC activities (\). \\ varied from 
0.03 to 1.0 (s _1 ). The same color curves dis- 
played in each panel were calculated from 
the same \ value. (B) The various kinetics of 
PI(4,5)P2 reduction (tFR) according to the 
changes in PLC activity (h) impact activation 
(left; A10-90%) and inactivation (right; A90- 
50%) of TRPC6/7 channel currents. Shown are 
the prediction made by the DG model (black) 
and the SPD model (red) for TRPC6 (circles) 
and TRPC7 channels (triangles). The star in 
the panels indicates the result from the h set- 
ting of 0.7. The solid and dashed lines indicate 
the current-rfK relationships experimentally 
observed in TRPC6 and TRPC7, respectively. 
The direction of the arrows and its size indicate 
the gap from the DG to SPD model. (C) Simu- 
lation results of PLGdependent characteristics 
of TRPC6/7 channels by DG (black) and SPD 
(red) models. (Left) The relationship between 
maximum open probability (P omax ) and Aj. 
(Middle) The time required to peak currents 
(s) and k. (Right) The total ionic influx and h- 
i=30 s. 
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activity increased in the range where the h values were 
small (k of <0.1). However, P omax in the SPD model was 
sustained or slighdy attenuated when the k values were 
>0.1. Kwon et al. (2007) showed that disruption of the 
potential PI(4,5)P2-sensing residue of the TRPC6 chan- 
nel nearly doubles the maximum current amplitude. 
This is consistent with our simulation, where the P om ax at 
k = 0.7 s eventually increased by two- to threefold in 
the DG model compared with that in the SPD model. 
The consistency of our simulation and the previous re- 
port strongly supports the functionality of the self-limiting 
regulation in TRPC6/7 channels for protecting the ex- 
cess cation influx. 

The second characteristic we measured was the shift 
in peak appearance. When we applied robust receptor 
stimulation, the time to peak currents were drastically 
shortened compared with the endogenous receptor 
stimulation (see PI(4,5)P2 dynamics at different levels... 
in Results) . Essentially, the similar tendency was repro- 
duced in both models (Fig. 8 C, middle) . However, the 
time to peak shortening was even more striking in the 
SPD model (Fig. 8 C, red symbols) than in the DG 
model (black symbols). Petersson et al. (2011) compu- 
tationally demonstrated that TRPC channel activity in 
the nervous system contributes to the temporal integra- 
tion of the generation of a lower rate-firing pattern at 
distal dendrites. Their simulation result suggests that 
the rate of firing may be controlled by the peak time of 
TRPC currents. Therefore, TRPC6/7 channels may play a 
critical role in determining the neural firing rate, de- 
pending on the strength of receptor stimulation. 

The third characteristic is the total ionic influx rela- 
tion (Fig. 8 C, right). The total ionic influx, which was 
calculated by integration of the flow during the recep- 
tor stimulation (t= 30 s), was predicted as the bell-shape 
response to the PLC activity in the SPD model. There- 
fore, the optimum PLC activity appeared neither weak 
nor robust, but was in the middle (0.1 < k < 0.3; Fig. 8 C, 
right, red symbols) . TRPC3/6/7 channels contribute to 
vascular tone, inflammation, cell remodeling, and intes- 
tinal contraction through increased Na + or Ca 2+ levels 
(Tsvilovskyy et al., 2009; Smedlund et al., 2010; Tauseef 
etal., 2012; Weissmann etal., 2012; Koenig et al., 2013). 
Hence, this bell-shaped relationship may be linked to 
these multiple pathophysiological responses. 

Collectively, channel inhibition by PI(4,5)P2 reduc- 
tion in DAG-sensitive TRPC3/6/7 channels not only ac- 
celerates the rate of current decay but also contributes 
to suppression of the current amplitude, shortening of 
the peak time, and the bell shape of the ionic influx, 
leading to the multiple functionality of these channels. 

PI(4,5)P 2 replenishment mechanism by diffusion 

It has been suggested that phosphatidylinositol 4-phos- 
phate is critical for replenishing PI(4,5)P2 via its phos- 
phorylation by PIP5K (Suh and Hille, 2002; Loew, 2007) . 



Furthermore, PIP5K substantially contributes to the quick 
recovery of PI (4,5) P 2 after its depletion by DrVSP activa- 
tion (Falkenburger et al., 2010; Imai et al., 2012). How- 
ever, in case of receptor-operated TRPC6/7 currents, 
we observed that a high concentration of ATP or its in- 
active analogue AMP-PNP in the patch-pipette solution 
had little effect on TRPC6/7 currents (Fig. S7), even in 
the plateau or biphasic phase. We speculate that an ad- 
ditional PI(4,5)P9 replenishment pathway, such as the 
detachment of PI(4,5)P2 from proteins, dispersion 
from its clustered complex (van den Bogaart et al., 
2011), translocation of phosphoinosi tides from the PIS 
organelle (Kim et al., 201 1) , or an unknown mechanism 
(Hammond et al., 2012), may be involved. 

However, here we focused on the lateral diffusion 
of PI(4,5)P2- When we analyzed FRET in regions of 
small compartments, a nonuniform reduction in FR 
by PI(4,5)P2 sensor upon receptor stimulation was ob- 
served (Fig. S3). This variability of a region-specific 
PI(4,5)P9 reduction was consistent with an electron 
microscopic study showing that PI(4,5)P2 located in 
a membrane microdomain (known as caveolae) de- 
creased slower than in an undifferentiated membrane 
upon receptor stimulation (Fujita et al., 2009). We 
have speculated that such nonuniformity of PI(4,5)P2 
reduction creates the opportunity for PI(4,5)P2 diffu- 
sion to occur. For that reason, we incorporated the 
pathway of PI (4,5) P2 replenishment by lateral diffusion 
into our model (Part 1 ) . By incorporation of this path- 
way, the plateau or biphasic currents after the peak cur- 
rent can be enhanced as the gap of PI(4,5)P2 dynamics 
between the local and global domain increases (Fig. S3) . 
In the case of TRP channels in flies, the plateau phase 
of the TRP current in response to light adaptation has 
been demonstrated to have a key role in phototransduc- 
tion (Lo and Pak, 1981; Minke, 2010). Therefore, incor- 
poration of replenishment of PI(4,5)P2 in our model 
may provide an important contribution to understand- 
ing the physiological consequences. 

Molecular insight into PI(4,5)P 2 and DAG binding 
to TRPC channels 

The fitting result of the SPD model predicted the dis- 
sociation constant of DAG binding to TRPC6 channel 
to be 20 to -40 pM and to TRPC7 to be <10 pM (Table 5, 
rows 17 and 18) . This parameter has not been well char- 
acterized, despite being critical for DAG-sensitive TRPC 
channel activation. In the experiment, receptor-operated 
TRPC7 currents demonstrated a higher sensitivity to DAG 
production (Fig. 3 E) than TRPC6 currents, which was 
consistent with the result predicted by the SPD model, 
thus providing a fundamental parameter into the TRPC6/ 7 
channel activation. However, it is still unclear exactly 
where these lipids bind to the TRPC3/6/7 channels. 

Among the TRP superfamily, heat- and capsaicin- 
activated TRPV channels are relatively well characterized 
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and have at least two PI(4,5)P2-binding sites in the 
C-terminal domains, proximal and distal regions to the 
S6 domain (Prescott and Julius, 2003; Doerner et al., 
2011; Ufret-Vincenty, et al., 2011). In the case of the 
TRPC family, earlier studies have also shown that the 
distal region of the TRPC6 C-terminal domain contrib- 
utes to the current decay (Kwon et al., 2007), whereas 
PI(4,5)P2 binding to the PH-like domain in the N-terminal 
region reduces TRPC3 current amplitudes (van Rossum 
et al., 2005). These results indicate that PI(4,5)P 2 bind- 
ing may be supported by multiple channel domains/ 
regions in a complex manner. Furthermore, another 
complex point regarding the PI(4,5)P2 regulation is 
that various heteromeric TRPC3/6/7 channels can be 
formed, including heteromers of TRPC3/6/7 (Goel 
etal., 2002) and other subfamilies (Hofmann etal., 2002; 
Takai et al., 2004; Ambudkar et al., 2006). Addressing 
these points is essential to further extend our under- 
standing on the self-limiting regulation by PI(4,5)P2 and 
DAG signals observed in TRPC3/6/7 channels. 
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